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Deriving minimum evolution times is of paramount importance in quantum mechanics. Bounds 
on the speed of evolution are given by the so called quantum speed limit (QSL). In this work we use 
quantum optimal control methods to study the QSL for driven many level systems which exhibit 
local two-level interactions in the form of avoided crossings (ACs). Remarkably, we find that optimal 
evolution times are proportionally smaller than those predicted by the well-known two-level case, 
even when the ACs are isolated. We show that the physical mechanism for such enhancement is due 
to non-trivial cooperative effects between the AC and other levels, which are dynamically induced 
by the shape of the optimized control field. 
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Introduction - The ability to precisely control quantum 
systems in a coherent way is one of the major challenges 
faced nowadays by physical sciences. This is because 
many technological developments, from electronic nan¬ 
odevices mm to quantum simulators and computers 
[21 m require unprecedented control over quantum 
systems. In this context, a lot of recent work has been 
devoted to the development of quantum control methods, 
both in theoretical mm and experimental IZIIHI grounds. 

An important part of theoretical quantum control 
analysis is the derivation of optimal evolution times 
[3110], mainly because it is desirable for quantum oper¬ 
ations to be implemented in the fastest possible way to 
avoid unwanted environmental effects which can destroy 
the coherence properties of the system. This task is 
usually tackled by quantum optimal control methods, 
either analitically [T0HI3] or numerically [HI [T5|. The 
typical problem in optimal control is to derive the field 
X(t) which optimizes a particular dynamical process 
(e.g. a transition from an initial state |'0o) to another 
goal state \tpg)) for a system described by a Hamiltonian 
H{X). 

The issue of time-optimal evolution is deeply con¬ 
nected with the fundamentals of quantum theory, as 
constraints on the speed of evolution are imposed by 
the time-energy uncertainty relation. The formal study 
of such constrains, due originally to Mandelstam and 
Tamm |16| . has led to the concept generally known 
as quantum speed limit (QSL) [TTHHj- The QSL 
determines the theoretical upper bound on the speed of 
evolution of a quantum system, and has been studied for 
both closed and open systems [23{2Sj- The connection 
between the time-dependent formulation of the QSL and 
quantum control has also been studied [57H23 ■ Although 
the usual formulation of the QSL can sometimes give no 
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meaningful bound for the evolution time in a quantum 
process, a recent work demonstrated a deep connection 
between the implementation of optimal control and 
QSL |30| . In particular, for a two-level system it was 
shown that the performance of the optimization is 
directly influenced by the relation between the (fixed) 
evolution time T fed to the algorithm and the quantum 
speed limit time Tqsl- the process did not converge if 
T < Tqsl- This result motivated the use of quantum 
optimal control as tool that allows to study the QSL 
in an heuristic fashion, specially for systems which are 
driven externally and are too complex to be studied 
analitically. 

In this letter, we implement quantum optimal control 
as a tool to study the QSL in many-level (i.e., more 
than two) systems, which show several avoided-crossings 
(ACs) in their energy spectrum. We show that the QSL 
time for processes involving more than one AC is in 
general smaller than the algebraic sum of the optimal 
times for each crossing, even when they appear to be 
well-isolated. We base our study on the analysis of 
a three-level system which shows two ACs. Through 
analysis of the obtained optimal fields, we study the 
physical mechanism involved in such speed up. We 
show that this field generates a non-trivial occupation 
of the levels which are initially uncoupled from the 
AC under consideration, but are on resonance with the 
applied field. Finally we extend our analysis to general 
A^-level systems, and show that the speed up reaches a 
saturation value as N increases. 

Three-level model and optimal control - We are 

interested in analyzing quantum systems which show 
ACs in their energy spectrum as an external parameter 
is varied. This model is of paramount importance in 
physics as it accounts for many interesting quantum 
phenomena, such as Landau-Zener transitions [3Ij . 
Landau-Zener-Stuckelberg interferometry [32] and quan¬ 
tum phase transitions [33] • Moreover, ACs can serve as 
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a pathway for designing robust control protocols [51H5S] . 
Real quantum systems usually present a complicated 
many-level spectrum with many ACs, and in most 
of the cases mentioned above, the analysis relies on 
the fact that one of such ACs can be regarded as isolated. 

We begin by considering the simplest extension of the 
two-level model showing one AC. Consider a three-level 
system described by the following Hamiltonian matrix 
given in the diabatic basis {|0), |1), |2)} 

/ 2A 0 \ 

HiX) = - 0 Ab \ = Ho + H(A) (1) 

^ V 0 Ab 2(A-eo) J 

where we define V (Hq) as the diagonal (non-diagonal) 
part of H{X). This model has been widely studied in 
quantum optics, as it is suitable for describing a three- 
level atom in a A configuration [37l[38]. In that context, 
the control parameters are usually A^ and Ab, which 
leads to the well-known STIRAP protocol for three-state 
control [3H]. Here, we take A to be our control parameter, 
which we allow to be time-dependent, while A^, Ab and 
Co are fixed parameters of the system. In Fig. (a) we 
show the energy spectrum of Hamiltonian Q as a func¬ 
tion of A. If eo ^ A^, Ab, the system presents two well- 
defined ACs, located at A = Aa = 0 (A = Ab = eg), with 
a mininum gap equal to A^ (Ab)- Under these condi¬ 
tions, it can be considered that only two states contribute 
to the dynamics of the system at each AC, while the re¬ 
maining one can be adiabatically eliminated m- Here, 
we will focus on control processes that connect some ini¬ 
tial diabatic state |^o) to a target The simplest 

example is depicted in Fig. [^(b), where |^o) = |1) and 
It/jg) = |0). As those states are involved in the AC la¬ 
beled “A”, it suffices to set the control field equal to Aa 
for a time = tt/Aa in order to generate the de¬ 
sired evolution. Moreover, state |2) does not take part 
in the process, and so this case could be described by a 
two-level system. Indeed, it has been shown m that, 
in this model, is indeed the minimum possible time 
required to connect two orthogonal states (as long as the 
minimum gap is fixed). We call the associated control 
field X{t) a “sudden-switch” field [35], which is a special 
case of the so-called composite-pulse protocol Uni EO], 
which is of bang-bang type (see also the discussion in 
[29|). We point out that this result can be interpreted 
in terms of the usual QSL formulation, which is derived 
from the time-independent Mandelstam-Tamm inequal¬ 
ity [l6] 


T> ^arccos{\{'ipo\i:{T))\), (2) 

where T is the evolution time AA^ is the variance of 
the Hamiltonian, AA^ = {{H — (iJ))^). An AC as 
discussed above can be described by a Hamiltonian 





Figure 1. (color online) (a) Eigenenergies of the Hamiltonian 
0 as a function of the control parameter A. The energy 
levels form two isolated avoided crossings at Aa and \b (for 
this example we used cq/Aa = 10 and As/Aa = 1). The 
eigenstates corresponding to each level (far from the ACs) 
are depicted with matching colours, (b) Schematics of control 
process I, which involves AC “A”. Here, the initial state ■ is 
|1) and the goal state © is |0). Below, we show a possible 
choice for the control field \{t) that generates such process, 
(c) Same as (b) but for process H, involving both ACs. The 
initial state ■ is |0) and the goal state © is |2). 


H = (A/2)cra, + Acr^, where ai are the Pauli matrices. 
Setting A = 0 for the sudden-switch field, and picking 
the initial state |'0o) as an eigenstate of dz, inequality 
( 0 ) gives precisely T > tt/A. 


A slightly more complicated process for the three-level 
system is depicted in Fig. [^(c). There, |V'o) = |0) and 
I'i/'g) = |2), which are not directly connected by any of the 
ACs in the spectrum. However, extending the previous 
protocol is straightforward: starting from |0), the con¬ 
trol field X(t) can be kept constant at Aa as before, and 
then switch instantly to Ab , where is kept constant for a 
time tt/Ab. In a such way, we can apply two succesive 
sudden-switch fields to generate the desired process. As 
a consequence, the total evolution time will be the sum 
of the required elapsed time at each individual AC, 



TT TT 

Aa Ab 


(3) 


Albeit being composed of two optimal protocols, it is 
not a priori clear that the multiple sudden-switch is also 
optimal. We recurr to optimal control methods in order 
to find the QSL time for this process. We use a Krotov 
optimization algorithm, as described in many previous 
works (see Refs. [IS] |13]). This procedure takes as an 
input the fixed evolution time T, an initial guess for the 
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Figure 2. (color online) (a) QSL time calculated from the 
optimal control procedure (see text for details) as a function of 
eo/Aa, for processes I (crossing one AC) and II (crossing two 
ACs). Dashed lines correspond to the time required by the 
sudden-switch protocol in each case, and . (b) Close- 
up of (a) in the region of low eq/Aa, where the ACs strongly 
interact, (c) QSL time for the case eq/Aa = 5 as a function 
of Ab/Aa, for process II. The dashed line corresponds to the 
prediction §. The same for the full line, but rescaled. 

control field A*^°^(f) and, of course, the Hamiltonian Q 
and both the initial and final states |^o) and \ipg). Each 
step k of the algorithm delivers an updated control field 
A(^)(t) and the corresponding evolved state (t)). 

We define the infidelity Ife = 1 — |^ 
as a measure of the success of the control process in 
the A:—th iteration. In our case, for every value of eq 
( the distance between the ACs), several runs of the 
algorithm were realized for different evolution times T. 
The initial guess was chosen in the following way. For 
each T we take the sudden switch field (see Fig. and 
appropiately shrink or expand it to fit each value of T. 
The field is then smoothed over to remove the discon¬ 
tinuities of the original function (the results we present 
are independent of the particular choice of smooth¬ 
ing), and also added an small correction linear in time. 
We will expand on the shape of the control fields later on. 

QSL time and optimal control fields - In order to 
obtain the desired QSL time for the model described 
above, we follow the method introduced in Ref. [30] . 
The basic idea is that the fidelity J^k = ^ cannot 
grow indefinitely if T should be smaller than the QSL 
time. Formally, for each value of T we look at the second 
derivative of (with respect to k) and analyze its sign. 
Then, the minimum value of T which gives 1''{k) < 0 
asymptotically, is chosen as the QSL time. 

We applied this procedure to the two control process 
discussed, and calculated Tq\b (* = 1, 2) in each case for 
several values of Eq, the parameter which measures the 
distance between the ACs in the energy spectrum (c.f. 


Fig. [Tl). The results we obtained are shown in Fig. ^ In 
the first case, it can be seen that is considerably 

larger than = tt/Aa for small eq, see Fig. H (b). 
This is, in principle, natural since in this regime the AC 
“A” is strongly affected by the AC “B”, which in turn 
leads to significant variations of the interaction rates 
(see Ref. m for more details). Away from that regime, 
TqgL converges to T^\ which is the well-known result 
for the two-level system, as discussed above. This result 
is indeed reasonable, since only the states |2) and |1) 
are involved in the process, and they both interact at 
AC “A”. However, it is interesting to point out that this 
behaviour allows us to quantitatively define the regime 
in which the ACs are well isolated. In the case shown 
in the figure, for which Ab = Aa, this is achieved for 
eo/Aa > 5. 

For process H, which involves both ACs, we first 
analyzed how the calculated QSL time changes when 
the size of one the gaps is modihed, for a fixed value of 
Eq. Results are shown in Fig. [^(c), where it can be seen 
that Tqg^ in fact scales as A^ , as expected from the 

dependance of the prediction Tg ' with the magnitude 
of the gaps, c.f. expression ([^. Indeed, the obtained 
data points can be fitted by Tggg = /3T^\ where 
/3 < I. This results implies that the calculated QSL 
time is smaller than Tg '. Remarkably, this result holds 
in all cases, even for large eq. The difference between 
TqgL and our prediction is larger for small eq, and 
decreases as the ACs are brought apart. However, for 
eo/Aa as large as 100, the difference is still larger than 
7%. This striking behaviour indicates that the control 
optimization can generate successful (i.e., with arbi¬ 
trary fidelity) control processes which are significantly 
shorter in time than the double sudden-switch, a pro¬ 
cess wich is time-optimal at each AC, as discussed above. 

In order to further analyze these results, we studied 
the optimal control fields derived from the quantum 
optimal control procedure, and the corresponding evolu¬ 
tions they generate. An example is shown in Fig. (a) 
and (b), for a particular value of eq/Aa- There, it can be 
seen that A(t) develops oscillations which are mounted 
on a step-wise function, the latter being a feature that 
is preserved from the initial guess. For all values of eq, 
the final shape of the field can be characterized mainly 
by the maximum amplitude A^ax and the frequency 
of the oscillations, which we obtain by calculating the 
fourier transform of A(t) . We plot those quantities as 
a function of the distance between the ACs in Fig. 

(c). Remarkably, we found that both scale linearly with 
Eq. Moreover, the equality fc = Eo/27r seems to hold 
perfectly for all cases considered. This analysis indicates 
that, in this model, the optimization leads to a control 
field that is not only reliable, but also simple and robust 
in the sense that it can be fully accounted for by means 
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^o/^B 


Figure 3. (color online) (a) Initial and optimized control 
fields X{t) for process II, using Ab/Aa = 1 and eo/ Aa = 10. 
(b) Time evolution of the populations Pk{t) for each one of 
the diabatic states \k} (k = 0,1,2), for the same parameter 
values as (a). Evolution time is fixed at T = Tq^j^ ~ 0.91Tg^^. 
The maximum value achieved by P 2 is shown, (c) Frequency 
fe (right axis) and maximum oscillation amplitude Amax (left 
axis) as a function of eo- Dashed lines indicate lineal depen¬ 
dences of both quantities with eo ■ 


of just a few paramaters, which show a well-defined 
dependance with the rest of the systems parameters, 
c.f. Fig. i (c). This result is indeed remarkable in at 
least two ways. In a pragmatic sense, it provides the 
means to develop a deeper theoretical understading of 
the time dependent problem. In fact, an explicit control 
protocol can be formulated and solved analytically for 
this problem [44] . But also, it provides further insight 
into the power of quantum optimal control theory. It 
is interesting to point out that, for a two-level model, 
optimal control techniques have been reported to lead 
(by extrapolation) to an analytical optimal control field 
m- However, for more complex systems, this procedure 
is expected to lead to complicated, non-analytical 
control fields. Also, we point out that some interesting 
discussions on the relation between complexity and 
quantum optimal control have been recently put forward 
in the literature IMlllT]. 

We will now describe how the particular shape of 
the optimized fields provides the physical mechanism 
for the overall speed-up of the evolution. As pointed 
out before, the field preserves the step-wise feature of 
the original sudden-switch field but with the novelty of 
showing oscillations. This behaviour can be understood 
in terms of navigation of the energy spectrum (Fig. [^. 


At the beginning of the protocol, the field activates 
the |0) e-?- |1) interaction by setting A = Aa. As a 
consequence, the initial population of state |0) starts to 
move to state |1). Simultaneously, the control parameter 
A starts to oscillate both rapidly (i.e. with a frequency 
/e ~ eo Aa) and with large amplitude (again, of 
the order of eg). This makes the control parameter A 
navigate areas of the spectrum close to the AC “B”. As 
a consequence, the previously uncoupled state |2) starts 
to get populated. Note that this is the target state \ijjg) 
for this process. When the first step is finalized, i.e. 
when the mean value of A switches to As, the three 
diabatic states have non-zero population. This is the key 
of the speed-up with respect to the sudden-switch field: 
the target state is already activated in the first half of 
the evolution. The protocol finishes by activating the 
dominant interaction |1) -o- |2), so that the population of 
12) continues to grow to 1. As happened in the first half 
of the evolution, the field also oscillates, now to drive A 
close to AC “A” with the purpose of depopulating state 
|1). It is interesting to point out that the simultaneous 
activation of multiple ACs seen in this process, can be 
regarded as an analogous effect to the wavefunction 
spreading reported in the optimal state transfer in spin 
chains [80] . 

Going through ACs in A-level systems - In the pre¬ 
vious paragraphs we showed how a process involving two 
ACs showed an enhancement of the QSL time. This re¬ 
sult motivates the question of wether this phenomenon 
can be further exploited by adding more levels to the sys¬ 
tem. In order to adress this problem, we generalize the 
three-level model described by the Hamiltonian ([T) to a 
A-level system. This can be done straightforwardly by 
defining the following Hamiltonian 


iV-1 
2 

Hn{X)= ^ (A-neo)|2n)(2n| 

n—0 

iV-2 1 

2 J 

+ ^ n eo|2n-I-l)(2n-I-1| 

n—0 

N-2 

+ + + + ( 4 ) 

n—0 

where [x] denotes the integer part of x. We schematically 
depict the energy spectra of these systems in Fig. [^ 
(a). Note that this models are constructed from the 
usual two-level spectrum with one avoided crossing by 
succesively adding diabatic levels which are parallel or 
perpendicular between each other. For A = 3, this 
scheme generates two crossings, but for A = 4 and 
A = 5, there are four and six crossings, respectively. 
We wish to study control processes involving (A — 1) 
ACs, as indicated in the Figure. To do so, we construct 
the Hamiltonian in such a way that only those crossings 
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Figure 4. (color online) (a) Schematic representation of 
some the A'^-level processes considered in this work. Gray 
lines correspond to the diabatic energy levels as a function of 
the control parameter. Blue lines show the path of the control 
process. The crossings between the levels which are traversed 
by the blue lines are avoided, while the others are exact (i.e. 
they do not interact, see main text for more details), (b) 
Ratio between the calculated QSL time Tqsl and as 

a function of N. Dashed line corresponds to the estimaton in 
ec. 0. The parameter r was chosen so to fit the data. For 
all cases, the distance between the ACs was set to eo = lOA, 
where A = 1 is minimum gap for all the ACs. 

are avoided, while the others must be kept closed (i.e. 
exact degeneracies in the spectrum). Note that there is 
no loss of genereality in making this assumption. If all 
the crossings are non-degenerate, the shortest control 
process which connects two-diabatic states will be the 
one which goes through the minimum possible amount 
of ACs (most likely, less than iV — 1). By closing some 
of the ACs, we are just forcing the system to follow a 
deliberate path. 

For this general model, we applied a procedure analo¬ 
gous to the one used for TV = 3 and extract Tqsl as a 
function of N. For simplicity, we show here results for 
the case in which any adjacent crossings are separated by 
a distance eo and all the gaps have equal sizes A„ = A for 
all n. Note that the sudden-switch protocol introduced 
at the beggining of this letter can be trivially adapted to 
this model. For a process going through N — 1 ACs in a 
A-level system, the total evolution time is proportional 
to N and given by 

Tf-'U(TV-l)^. (5) 

The hgure of merit for this analysis can be defined as 
the ratio between the calculated QSL time and 
We plot this quantity for several values of N in Fig. 
(b), where it can be seen that it decreases as N grows. 


This remarkable result shows that the control processes 
benefits from the existence of more levels. It is even more 
interesting to note that this result can be derived just by 
considering the speed-up already registered in the N = 3 
case. As seen above, the key to the speed-up is that in 
each step, the control field draws a small population to 
an adjacent level, originally not coupled to the AC. Also, 
given the shape of the spectrum, there cannot be more 
than three ACs simultanously active during the protocol: 
the one activated by the sudden-switch field and its ad¬ 
jacent neighbors. As a consequence, the QSL time will 
be smaller or equal than for all A > 2, and the 

difference is proportional to N, so that we can write 

Tqsl = - (A - 2)r (6) 

for some value of r < Tg^ This leads to 


The parameter /3(A) is then found to decay with A 
to a constant value below 1, thus proving the existence 
of a speed-up for this control process. 

Final remarks - In this letter we studied the QSL of 
systems which show more than one (localized) AC in 
their energy spectrum. For this purpose we exploited 
quantum optimal control methods to find the minimum 
evolution time required to implement control processes 
involving many ACs. We found that the procedure 
introduced in Ref. |30j allows to correctly discriminate 
the QSL time for A-level systems with multiple ACs. By 
analizing two different control processes in a three-level 
model, we show the minimum time required to cross 
two ACs is smaller than the sum of the optimal times at 
each crossing. We then generalized this result to A-level 
systems, showing that the derived speed-up still holds 
but reaches a saturation points as A increases. As an 
interesting byproduct of our analysis, we found that 
the optimized control fields derived in this procedure 
are easily described by means of a few well-defined 
parameters. Thanks to that, we identihed that the 
physical mechanism that cause the QSL enhancement is 
based on the collective effects between the states which 
interact at each AC and others, which are dynamically 
coupled by the control fields. For future studies, it would 
be interesting to relate the observed phenomenon to the 
recently demonstrated environment assisted speed-up of 
the evolution of a quantum system |45| . In this context, 
the systems initial |V'o) and final states \ipg) could be 
interpreted as forming a two-level system interacting 
with an environment which produces the aforementioned 
speed-up via a non-markovian effective dynamics. 
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